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We describe a mode sorter for two-dimensional parity of transverse spatial states of light based on 
an out-of-plane Sagnac interferometer. Both Hermite-Gauss (HG) and Laguerre- Gauss (LG) modes 
can be guided into one of two output ports according to the two-dimensional parity of the mode in 
question. Our interferometer sorts HG n m input modes depending upon whether they have even or 
odd order n + ra; it equivalent ly sorts LG l p modes depending upon whether they have an even or odd 
value of their orbital angular momentum /. It functions efficiently at the single-photon level, and 
therefore can be used to sort single-photon states. Due to the inherent phase stability of this type 
of interferometer as compared to those of the Mach-Zehnder type, it provides a promising tool for 
the manipulation and filtering of higher order transverse spatial modes for the purposes of quantum 
information processing. For example, several similar Sagnacs cascaded together may allow, for the 
first time, a stable measurement of the orbital angular momentum of a true single-photon state. 
Furthermore, as an alternative to well-known holographic techniques, one can use the Sagnac in 
conjunction with a multi-mode fiber as a spatial mode filter, which can be used to produce spatial- 
mode entangled Bell states and heralded single photons in arbitrary first-order (n + m = 1) spatial 
states, covering the entire Poincare sphere of first-order transverse modes. 



I. INTRODUCTION 

The manipulation of the transverse spatial degrees of 
freedom of photons offers avenues for quantum infor- 
mation studies using single photons as qubits. Several 
recent studies have focused on higher order transverse 
spatial modes, of which the Hermite-Gauss (HG) and 
Laguerre- Gauss (LG) modes are examples In particu- 
lar, several different interferometers of the Mach-Zehnder 
type have been experimentally employed for sorting sin- 
gle photons by their transverse mode properties. For the 
LG modes, an orbital angular momentum (OAM) sorter 
was demonstrated at the single-photon level which sorted 
modes with even OAM values, which exited one port of 
the interferometer, from those with odd OAM values, 
which exited the other port @, @]- For the HG modes, a 
mode index sorter was theoretically proposed to measure 
the value of a single HG index norm [J] , while a different 
setup experimentally distinguished the HGio and HGoi 
modes, thereby acting as a transverse-mode beam split- 
ter for first-order HG modes, in analogy with a polarizing 
beam splitter A similar setup was later used to per- 
form various manipulations on photonic qubits using a 
one-dimensional transverse spatial degree of freedom [6[ . 

The two experimentally implemented sorter designs 
and H, shown in Fig. [TJ share several similarities. 
Among these are their Mach-Zehnder character with two 
independent arms, and their reliance on complete con- 
structive and destructive interference at the output beam 
splitter in order to sort certain transverse spatial modes 
into one of two distinct output ports. Also, they both 
have the advantage of not relying on computer-generated 
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FIG. 1: Previously realized 1-D and 2-D parity sorters of the 
Mach-Zehnder type, each with a tilt able phase-shifting glass 
plate (GP) in one arm. (a) The 1-D parity sorter is distin- 
guished by having an extra mirror in one arm. Upon enter- 
ing this interferometer, an LG oi mode is sorted into its con- 
stituent HGio and HGoi components due to complete con- 
structive and destructive interference at the output ports, (b) 
In contrast, the 2-D parity sorter has a Dove prism in each 
arm. One of the Dove prisms is rotated 90° with respect to the 
other, which causes a 180° relative rotation of the two inter- 
fering beams in the transverse plane. In this case, both HGio 
and HGoi modes exit the same port, so that an incident LGo\ 
mode is not decomposed. 



holograms, which are lossy in practice, so that the sort- 
ing efficiency was limited only by the transmission and 
reflection coefficients of the constituent mirrors, glass 
plates, and dove prisms. However, though they are simi- 
lar, the two aforementioned sorters actually treat trans- 
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FIG. 2: A tabulation of the action of the the 1-D sorter upon 
the HG modes of order n + m < 3. The plus signs designate 
modal lobes which are in phase with each other and 180° out 
of phase with the unmarked lobes. Modes with an even value 
for n exit port A, while modes with odd n- value exit port B 
(see Fig. |H (a)). 



FIG. 3: A tabulation of the action of the the 2-D sorter upon 
the HG modes of order n + m < 3. The plus signs designate 
the relative phases of the lobes as in Fig. [2] Modes with 
an even value for n + m exit port A, while modes with odd 
n + m- value exit port B (see Fig. [T] (b)). 



verse spatial modes in fundamentally different ways. The 
first order HG mode beam splitter in Fig. [T] (a) distin- 
guishes modes on the basis of one-dimensional (1-D) par- 
ity about the axis normal to the plane of propagation, due 
to the extra mirror reflection in one of the interferometer 
arms. This means in particular that HGio and HGoi 
input modes propagating in the x-z plane exit different 
ports as shown since they respectively have odd and even 
parity with respect to reflections about the ?/-axis, de- 
scribed by the transformation E (x, y) — > E (— y). On 
the other hand, the OAM sorter in Fig. [T] (b) distin- 
guishes modes on the basis of the ^o-dimensional (2- 
D) parity of the transverse (x-y) plane-that is, whether 
the mode is even or odd with respect to the transfor- 
mation E (x, y) — > E (—x, —y). Therefore, conversely to 
the prior case, the HGio and HGoi input modes will be 
sorted into the same port by the OAM sorter. Since, 
to our knowledge, no experiments have been reported in 
which a sorter of the OAM type was employed to sort 
HG modes, this fact seems to have been overlooked and 
the difference between the two sorters has not been em- 
phasized until now. We will henceforth refer to sorters of 
the aforementioned types as 1-D and 2-D parity sorters, 
respectively. 

Figs. [2] and [3] give a general outline of how HG nrn 
modes are sorted by the 1-D and 2-D parity sorters. Since 
the reflection from the extra mirror of the 1-D sorter 
inverts the sign of the x-variable only, it sorts according 
to whether the single index n (which corresponds to the x 
spatial variable) is even (output port A) or odd (output 
port B). For the 2-D sorter, the two interfering beams 
undergo a 180° relative rotation in the transverse plane 
(equivalent to sign inversion about both the x and y axes) 



so that it sorts according to whether the sum of the two 
indices n + m is even (port A) or odd (port B). 

A major disadvantage of both the 1-D and 2-D parity 
sorters mentioned above is their Mach-Zehnder charac- 
ter, or more specifically the existence of two indepen- 
dent beam paths in each interferometer. This makes 
them subject to phase noise and drift so that in prac- 
tice it is difficult to keep them aligned. In this work, we 
overcome this difficulty for the latter design by demon- 
strating experimentally a novel type of 2-D parity sorter 
in the form of an out-of-plane Sagnac interferometer. 
Since the two interfering paths of the Sagnac are counter- 
propagating (see Fig. the interferometer is automat- 
ically phase-stable. As a result, experiments involving 
multiple cascading interferometers are considerably more 
feasible when employing a Sagnac then with a Mach- 
Zehnder interferometer. Such cascading allows full 2-D 
sorting of arbitrary superpositions of modes, including 
the OAM sorting schemes described in 0, @1 • 

In what follows we discuss the design of our out-of- 
plane Sagnac interferometer. We present, to our knowl- 
edge, the first direct 2-D parity measurements of HG 
modes. We also propose several applications of this in- 
terferometer, including its use as an alternative to holo- 
grams in spatial mode filtering, the production of Bell 
states entangled in first-order transverse spatial modes, 
and the production of heralded single photons in first- 
order transverse spatial states corresponding to an ar- 
bitrary point of the first-order spatial mode Poincare 
sphere [7[. 
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FIG. 4: The 2-D parity sorting Sagnac interferometer. Mir- 
rors M1-M3 cause the beam path to travel out of the x-z 
plane, tracing out two congruent sides of an isosceles triangle 
in the x-y plane, with congruent base angles and apex angle 
P — TV — 20. Incident HGio and HGoi modes exit the same 
port of this interferometer, as do LGoi modes, as is the case 
for the sorter in Fig. [T](b). See text for further discussion. 



II. PHASE-STABLE PARITY SORTING 
SAGNAC INTERFEROMETER 

The Sagnac interferometer can be made to operate as 
a mode sorter because of a geometric phase accumula- 
tion implemented by the out-of-plane reflections at three 
of the interferometer mirrors, which results in a rotation 
of the transverse spatial profile of the beam through an 
angle ft about the beam axis [8, 9]. To see this, consider 
a simplified version of our 2-D parity-sorting interferom- 
eter, shown in Fig. [H The mirrors labeled M1-M3 in 
the figure are oriented such that the reflected beam path 
travels out of the x-z plane, and traces out the two con- 
gruent sides of an isosceles triangle which lies completely 
in the x-y plane. The action of the 50:50 beam splitter 
(BS in the figure) gives rise to two counter-propagating 
paths of equal path length, which is the reason for the 
phase-stability of the interferometer. We denote the path 
involving reflections off of the mirrors M1-M4 in the suc- 
cessive order M4-M3-M2-M1 as the clockwise path, and 
that involving reflections of order M1-M2-M3-M4 as the 
ant i- clockwise path. 

A convenient parameter for discussing the state of the 
beam is the unit helicity vector h, defined as h = k 
where n is the number of reflections, and k is the beam 
wave vector Q. Employing this, we now consider the four 
vectors hi through I14 associated with the anti-clockwise- 
traveling beam, defined as the unit helicity vectors prior 
to the respective reflections off of M1-M4 as shown in 
the figure. These vectors define four points P1-P4 on 
the unit sphere (See Fig. [5] (a)); however, due to the 
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FIG. 5: (a) The unit sphere construction corresponding to 
the beam path in Fig. 2] as discussed in the text. The beam 
Dath angles and (3 from Fig. 2] define the helicity vectors 
hi — I14, which in turn define the points P1-P4 on the spheri- 
cal surface as shown. The dashed lines denote the triangular 
closed loop that results from successively connecting these 
points with geodesic curves. The figure is drawn for the spe- 
cial case where = 45°, so that /3 = 90°. In this case, the 
angles between the distinct helicity vectors are all 90°, so that 
the area enclosed by the loop is exactly one eighth of the area 
of the entire sphere, which corresponds to a transverse rota- 
tion of O = ^ rads, or 90°. (b) Plot of the beam rotation 
angle Q vs. the parameter 0. (c) Plot of the relative rotation 
angle ^ (modulo tt) of the two counter-propagating beams 
after passing through the Sagnac interferometer vs. 9, em- 
phasizing the two special cases where ^ = 180° and ^ = 90°. 
These two cases correspond to the first and second cascaded 
interferometer stages for the OAM sorting scheme in Q. 

odd number of out-of-plane reflections, we have hi = 
£14, which implies that PI = P4, so that only three of 
these points are unique. Therefore, connecting in order 
the points P1-P4 on the spherical surface with geodesic 
curves (great circles) results in a triangular closed loop 
as shown in the figure. 

The spherical surface area enclosed by this loop is equal 
to the azimuthal rotation angle Q of the transverse beam, 
as discussed in [8[. In the present special case of three 
mirror reflections resulting in a closed triangular geodesic 
surface, one may apply the following relation (due to Eu- 
ler) involving the surface area ft of the sphere [TO]: 

Q 1 + cos a + cos (3 + cos 7 

cos — = -5 (1) 

2 4 cos t| cos ^ cos ^ 

In (pQ), a = hi • ti2, (3 = I12 • h3, and 7 = I13 • are the 

angles between the respective helicity unit vector pairs. 
In our present case (see Figs. H]|5]), a = 7 = ^, and 
/3 = tt — 20, so that the right hand side of (pQ) simplifies 
to 1+cos j> 7 ^~ 20 ? , which is equivalent to sin#. Therefore, 

2cos(f-0) ' 1 ' 



we find for our design that Q depends on through the 
simple relationship 

cos — = sin# (2) 

The plot of Q vs. in Fig. 03(b), which has been com- 
posed from ([2j), shows that Q can be tuned to any angle 
between 0° to 180° by choosing the appropriate value of 
the parameter which ranges from < < it. Now, 
for the clockwise traveling beam, the direction of the 
geodesic path of Fig. [5ja) is reversed, so that the trans- 
verse rotation is through angle — O, or opposite that of 
the anti-clockwise traveling beam. We therefore find that 
the absolute relative rotation \I/ between the two counter- 
propagating beams can be expressed as ^ = tt— |2fi — 7r|, 
as given through the plot of ^ vs. in Fig. 0(c). 

By orienting Ml and M3 in order to select the succes- 
sive values = j , ^ , ^ , . . . one can achieve the relative 
rotations ^ = 7r, ^, J, ... required for each of the inter- 
ferometer stages of the OAM sorter scheme in 0, E| . As 
discussed in [3], a tiltable phase-shifting glass plate in 
one of the Mach-Zehnder arms introduces an adjustable 
relative phase difference between the two interferometer 
paths, in addition to the phase difference already intro- 
duced by the transverse spatial rotation. This additional 
relative phase is needed in order to distinguish between 
modes possessing odd OAM values without the use of 
holograms, which are lossy in practice. Since our device 
is a common-path interferometer, a glass plate cannot 
introduce the required relative phase shift between the 
counter-propagating interferometer paths. However, one 
may surmount this difficulty by employing a device con- 
sisting of a tiltable birefringent waveplate surrounded on 
both sides by Faraday glass which produces the requisite 
phase shift as discussed in Appendix [A[ Therefore, like 
the device of Fig. [T](b), our interferometer is able to sort 
and measure photons possessing arbitrary absolute OAM 
values with 100 % efficiency in principle. Furthermore, as 
our Sagnac interferometer is phase- stable, it provides, for 
the first time, a practical way to realize a multiple-stage 
cascaded sorting scheme for single photon states. 

In the special case where 6 = 45° is chosen, the rel- 
ative rotation ^ = 180° and our Sagnac interferometer 
becomes a 2-D Parity sorter. We note that this type of 
Sagnac interferometer (with = 45°) was used in Q to 
measure the transverse spatial Wigner function for light 
at the single-photon level, while a similar interferome- 
ter was previously employed for similar purposes in [ill ]. 
The fact that the transverse Wigner function is the ex- 
pectation value of the 2-D parity operator displaced in 
phase space [q|, [l2| provides another way to look at our 
sorter, which simply measures the non-displaced parity 
operator (corresponding to the transverse Wigner func- 
tion evaluated at zero transverse position and momen- 
tum) for the special case of 6 = 45°. Following [11], one 
could in principle use a single rotatable Dove prism in the 
place of mirrors M1-M3, however there are various disad- 
vantages to this approach including a significantly lower 
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FIG. 6: The Experimental Apparatus. The propagating and 
counter-propagating beams interfere at a 50:50 beam split- 
ter (BS1), while a second 50:50 beam splitter (BS2) sepa- 
rates the backward-propagating output mode (port A) from 
the forward-propagating input mode. An external cavity he- 
lium neon laser (HeNe) acts as the source, with thin crossed 
wires (W) and an iris (I) inserted into the cavity in order 
to select the higher order HG modes. A third beam splitter 
(BS3) picks off part of the source beam to use a a reference 
beam (Ref) for the interferometry experiments discussed in 
the main text. Two Berek polarization compensators (BC) 
are placed inside the interferometer to correct for the Fresnel 
polarization changes due to reflections from the out-of-plane 
dielectric mirrors. The output fields at both port A and port 
B were measured with CCD cameras. 



numerical aperture and the introduction of astigmatism 
[9]. One may avoid these difficulties by substituting the 
out-of-plane mirror configuration described in [13]. Fur- 
thermore, since an "all mirror" Sagnac has already been 
successfully employed (for other purposes than the sort- 
ing of HG or LG modes) at the single-photon level in 
[9[, we conclude that a Sagnac with = 45° is capable 
of measuring the 2-D parity of single photon states (or 
equivalently, whether a single photon's OAM is even or 
odd). 

Finally, we note that our 2-D parity sorter can also 
be used to sort the output of a three-mode fiber into its 
constituent zero- and first-order mode families with 100% 
efficiency in principle. In addition to presenting direct 2- 
D parity measurements, we carry out this fiber filtering 
experiment in Section [Till an d discuss applications of this 
technique to quantum information processing in Section 

ED 



III. 2-D PARITY SORTING EXPERIMENT 

Our experimental apparatus was set up as shown in 
Fig. [6l which includes two Berek polarization compen- 
sators as well as two additional 50:50 beam splitters in 



Experiment 



Theory 



Input 

Mode 



HGo 



HGo 



HGi 



Input 



Interference 



Intensity Port A Port B Intensity 



Profile 



Profile 



HG 4 



HGu 



HGi 



HG 3 



3 Mode 
Fiber 




Input 

Mode 



HGo 



HGo 



HGi 



HG45 



HGi 



HGi 



HG 3 



3 Mode 
Fiber 



Input 
Intensity 
Profile 



Port A 




1 


1 


II 


II 


t 


• 



nan 



FIG. 7: Observed output intensity profiles and interference 
patterns for a given input field. 



FIG. 8: Predicted output intensity profiles and interference 
patterns for a given input field. 



relation to the simplified version of the Sagnac shown in 
Fig. [U The Berek compensators (New Focus model num- 
ber 5540) were inserted into the interferometer as shown 
in the figure in order to improve the interference visi- 
bility by correcting for the Fresnel polarization changes 
due to the reflections from the dielectric mirrors M1-M3. 
Without the compensators, the propagating and counter- 
propagating beams ended up with unequal elliptical po- 
larization states upon interference at the beam splitter, 
which led to poor interference visibility and therefore in- 
efficient sorting. Since one of the interferometer's output 
ports (port A in Figs. [H and [6]) is counter-propagating 
with respect to the input beam, the addition of beam 
splitter BS2 acted to deflect half of the output beam from 
this path so that it could be imaged on the CCD as shown 
in the figure. 

Thin crossed wires were inserted into the laser cavity 
in order to excite higher order HG modes by suppress- 
ing lower order transverse modes which lacked a nodal 
point near the location of the wire. An intracavity iris 
was then constricted in order to suppress the remaining 
excited modes that had a higher order then the desired 
mode. For a given input mode, we measured the out- 



put intensity profiles at both port A and port B of the 
interferometer using a CCD camera. In order to demon- 
strate the phase structure of the sorted HG modes, we 
employed a third beam splitter (BS3 in the figure) in or- 
der to pick off a portion of the input beam to use as a 
reference beam. We then interfered the reference beam 
for several input modes with their corresponding parity- 
sorted output modes by steering the input and output 
beam paths together so that they were co-propagating 
with a slight misalignment, and recorded the resulting 
intensity patterns. In a second experiment, the wire and 
iris were removed and the light was coupled into and 
out of a three-mode fiber (Thorlabs SM-780) which then 
acted as the input source. 

We note here that due to the 50% loss of the signal 
from port A from the presence of the beam splitter BS2, 
this setup was not capable of sorting single photons with 
high efficiency. However, by replacing beam splitter BS2 
with a Faraday isolator and appropriately adjusting the 
compensators, 100% efficiency can be obtained in princi- 
ple as shown in Appendix iBl 

A representative sample of our results are presented 
in Fig. [7J which agree with the corresponding theoreti- 



cal predictions in Fig. [HJ In these figures, the first col- 
umn indicates the input mode, while the second gives the 
observed/predicted intensity distribution for that mode. 
Columns three and four give the observed/predicted dis- 
tributions at output ports A and B, respectively. The 
final column gives the observed/predicted interference 
patterns between the input and output modes as dis- 
cussed above, which exhibit a characteristic "forking" of 
the vertical interference fringes which shows the phase 
structure of the modes [23[. Note that the the HGqq 
and HGn modes exit port A, while the HGio an d HGqi 
modes exit port B regardless of their orientation as pre- 
dicted in Figs. [3] and [U Furthermore, the HG15 mode 
(which has even parity since n + m = 1 + 5 = 6) ex- 
its port A while the HG32 mode (which has odd parity 
since n + m = 3 + 2 = 5) exits port B. We note here 
for completeness that a p th order LG l ±p mode is sorted 
by this interferometer according to whether p is even or 
odd which parallels the current demonstration of HG nm 
mode sorting according to the even/odd parity of n + m. 
We have therefore shown with these results that stable, 
cascadable, single photon OAM sorting schemes corre- 
sponding to those discussed in [2|, [3[ are indeed feasible 
with our interferometer. 

Finally, we note that the bottom row of Fig. [71 shows 
that the output field from the three-mode fiber is de- 
composed into the fundamental HGqq mode, which exits 
out of port A, and an equal superposition of first or- 
der Laguerre-Gauss fiber eigenmodes, which interfere to- 
gether to make a first order HG mode which exits out of 
port B (the theory plot in Fig. [8] is for a field composed 
of 85% HGoo and 15% first order modes). This explic- 
itly demonstrates the usefulness of our 2-D parity sorter 
as a "beam splitter" for even and odd transverse spatial 
modes, while at the same time suggests that it can be 
employed as an alternative to relatively lossy holograms 
in spatial mode filtering, as discussed below. 



IV. APPLICATIONS 

A. First Order Mode Filter 

One particularly useful characteristic of our Sagnac in- 
terferometer is its ability to distinguish between zero- 
order and first-order transverse HG modes. When used 
in conjunction with a three-mode optical fiber, this abil- 
ity allows the Sagnac to act as a spatial mode filter that 
passes only the first-order HG modes while rejecting all 
other mode orders. To see this, consider a linearly po- 
larized monochromatic paraxial beam with electric field 
E (x, y, z) e~ ZUJpt coupled into an optical fiber with input 
and output coupling lenses as shown in Fig. [9] (a) (we 
align the fiber with the z-axis and place the origin at 
the fiber input face). Such a beam comprises a general 
solution to the scalar paraxial wave equation, of which 
the paraxial HG modes form a complete basis [14]. Em- 
ploying the notation E (x, y,z) = \^o), one may therefore 
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FIG. 9: A spatial filter for zero and first-order HG modes, 
(a) A three-mode optical fiber (3MF) with input and out- 
put coupling lenses (CL). The fiber, which has a parabolic 
refractive index profile, acts as a mode filter upon an ar- 
bitrary transverse input beam |^ ) = X^ m =o Cnm \HG nm ) 
such that upon coupling into and out of the fiber the resulting 
beam state is of the form \i/jf) = coo \HG00) + cio \HG 10) + 
coi \HG01). (b) Inserting a 2-D Parity Sorter (S) after the 
fiber separates the zero and first-order modes into ports A 
and B, respectively. 



write 

00 

|V>o> = E Cnm \ HG ™) ( 3 ) 

n,m=0 

where c nm = (^0 \HG nm ) is the overlap integral between 
the incident field and the HG nm mode, and the constant 
Eq is the electric field magnitude. When the fiber and 
input coupling lens are properly aligned, the beam will 
be centered coaxially with the fiber, and the location of 
the beam waist will coincide with the plane of the fiber 
input face at z = 0. Under these conditions, the paraxial 
HG modes evaluated at the fiber input plane have the 
simplified transverse form 

\HG nrn ) 

= A nm — H n (V2—) Em (V2— ) e~^~^ (4) 



where A nm = y 2 ™+"Wm! > H p ( u ) in ® denotes a Her- 
mite polynomial of order p, and wq is the beam waist 
radius (([4]) is obtained by taking the limit as z tends to- 
ward zero of the associated paraxial expression of [141]). 

In comparison, the eigenmodes of an optical fiber with 
a parabolic refractive index profile n(p) = n (1 - A^J 
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obeying the weakly guided condition A <C are well ap- 
proximated |l5| by the linearly polarized Laguerre-Gauss 
modes \LG q ^ m ): 



\LG. 



i,m) - ( P a ) 



\m\ 



L \rn\ 



v p - 

a 2 



VP ■ ± 



(5) 



where a is the fiber radius, V 

-\m\ 



/ zA is the normal- 



ized frequency, the L\ are generalized Laguerre poly- 
nomials, and we have employed transverse cylindrical 
coordinates (p, 0). Furthermore, for sufficiently small 
a, the higher order modes will be cut off [16[ so that 
the fiber will support only the three LG modes |LGo,o)> 
|LGo,+i), and |LGo,-i), or equivalently the three HG 
modes \HGoo), \HGiq), and \HGq\). Therefore, regard- 
less of the form of the transverse input state |^o) in equa- 
tion (|3|) , a weakly guided parabolic fiber will in principle 
filter out the higher-order HG mode contributions while 
transmitting the fundamental and first-order modes with 
high efficiency such that the remaining output state ipp 
will have the form 
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FIG. 10: Possible applications of a 2-D parity sorter to quan- 
tum information processing, (a) Proposed scheme to produce 
Bell states entangled in first-order transverse spatial modes, 
(b) Proposed scheme to produce heralded single photons in 
arbitrary first-order transverse spatial states. 



\^f) = coo \HGoo) + cio \HG 10 ) + c 01 \HG 01 ) (6) 



provided that the mode matching condition wq — a 




is met. Given this fiber-filtered state, our Sagnac can 
then sort \^f) such that the output modes at port A and 
port B are \ip A ) = c 00 \HG 00 ) and \^ B ) = cio \HGio) + 
coi |iJGoi), respectively, as shown in Fig. [9](b). 

B. Applications to Quantum Information 
Processing 

Here we propose two applications of 2-D parity sorters 
to quantum information processing (QIP) with trans- 
verse spatial modes. Both make use of spontaneous 
parametric down conversion in addition to the first-order 
mode filter discussed above, as shown in Figs. flOl(a) and 
(b). In both cases, a crystal with a nonlinearity is 
cut and oriented for Type-II collinear phase matching and 
pumped with a sufficiently weak cw laser in a well-defined 
HG mode HG nmi such that the resulting state is predom- 
inantly a superposition of the vacuum and a two-photon 
state |^) = Ci 1 1 vac)) + G 2 \\^2 P h)), where C\ > G 2 and 
the "double-ket" symbol || )) denotes fock-space states 
of the electromagnetic field [17[ . 

When sufficiently narrowband filters centered at half 
the pump frequency uo p are employed, the co-propagating 
two-photon state \^2ph) is frequency degenerate, and is 
therefore not entangled in the spectral degree of free- 
dom. However, photons have four degrees of freedom, 
and Y&2ph) is entangled in each of the remaining three 
(polarization and the two transverse coordinates). In 
particular, if the appropriate birefringent phase compen- 
sators are used [18|, then the biphoton mode function 



corresponding to the two-photon component \^2 P h) of 
the state |^) can be written as pH [T^ 

\*hph) = \4>nm) ® ( \H) \V) + \V) \H) ) (7) 

where \H) and \V) denote the horizontal and vertical 
single-photon polarizations, and \(j) nm ) denotes the spa- 
tial portion of the biphoton mode function. In [19[, \4> nm ) 
was expanded in terms of the HG modes for an arbitrary 
HG pump beam HG nrni and shown to be entangled: In 
terms of the HG mode functions, |0 nm ) takes the form 

oo 

\4>nm)= E c^ st \HG jk )\HG st ) (8) 

j,k,s,t=0 

where the specific coefficients CjJ^ t depend on the order 
nm of the pump beam via equation (18) of [19]. As 
shown there, the expansion contains an infinite number of 
nonzero terms, and converges slowly in general. However, 
many of the coefficients c™™ st are zero due to selection 
rules. We apply this result below in two special cases. 

1. Bell State Generation 

The first QIP application we discuss is the production 
of Bell states entangled in first-order transverse spatial 
modes, as shown in Fig. [10] (a). An HG Bell state gen- 
eration experiment similar to this has been carried out 
by [10], which used holograms to sort and detect spatial 
modes. However, employing our Sagnac interferometer 
in the place of holograms will provide an improvement in 
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efficiency over standard holographic techniques. In the 
current setup, the pump beam is in an ifGoo mode, and 
following [19] we find that the spatial part of the down 
converted biphoton mode function has the form 

\M = c \HG 00 ) \HGoo) 

+ ci ( \HG W ) \HG W ) + \HG 01 ) \HG 01 ) ) 

+ c 2 ( | if G 00 ) \HG 02 ) + \HG 02 ) \HGoo) 

+ \HG 00 ) \HG 20 ) + \HG 20 ) \HG 00 )) 

+ ... (9) 

where c\ = |co ~ 0.04, and C2 ~ —0.03 under typical 
experimental conditions (pump beam width 0.1 mm and 
crystal length 1 mm), and all other coefficients of total 
order j + k + s + t = 2 are zero. The omitted terms 
in ([9]) all involve biphoton mode functions of total order 
j + fc + s + t > 4. However, since each individual photonic 
mode is filtered by the Fiber-Sagnac combination as dis- 
cussed above, the mode function in (0 retains only its 
first-order component, so that the (renormalized) mode 
function exiting port B is given by © and ([9]) as 

|Vb> = ^= ( \HG W ) \HG W ) + \HG i) \HG i)) 

®(\H)\V) + \V)\H)) 

(10) 

as shown in the figure. Thus, the filtered biphoton mode 
function is maximally entangled in both polarization and 
first-order spatial modes. This is a Bell state. 

We note here that the two photons comprising the Bell 
state are co-propagating; in order to turn this source into 
"useful" entanglement, one must separate these photons 
onto two separate paths without destroying the entangle- 
ment relationship. This can be done by inserting a 50:50 
beam splitter after port B, albeit with a 50% loss of the 
entanglement source. To obtain the HG-entangled Bell 
states without such a loss, one may instead use a polariz- 
ing beam splitter. The co-propagating photons are then 
path-separated with the cost of destroying the polariza- 
tion entanglement, while the spatial mode entanglement 
remains. 



2. Heralded Photons in Arbitrary Spatial Modes 

We now consider the second application of our mode 
sorter, shown in Fig. [10] (b). In this setup, the pump 
beam is a first-order HG\q mode rotated through an an- 
gle of 45° in the transverse plane, denoted as HG450. 
Since HG450 comprises an in-phase superposition of 
HG10 an d HGqi modes, we express it in terms of these, 
and again employ equation (19) of [19[, finding that that 
down- converted biphoton mode function is 



1045°) = |0io) + |0Ol) 

= ci ( \HG 10 ) \HG00) + |# Goo) \HG 10 ) 

+ \HG01) \HG00) + |# Goo) \HG i) ) + ... 

(11) 

where all other coefficients of total order j + k + s-\-t = 
1 are zero. Similarly to the previous case, the omitted 
terms in ([TTj) all involve biphoton mode functions of total 
order j + k + s + 1 > 3, and are therefore filtered by the 
fiber. Therefore, the action of the Sagnac on the present 
fiber-filtered biphoton mode function is to path-separate 
the photons in the fundamental mode 1 00) from those in 
the higher modes. Thus, upon exiting the Sagnac the 
biphoton mode function is 

M = \HG o) A |#G 45 o) B ( \H) A \V) B + \V) A \H) B ) 

(12) 

where \HG450) = \HG 10 )^\HG 01 ) and A and B label the 
output ports, so that the detection of a | if Goo) photon 
in port A heralds a single polarization-entangled photon 
in the pure spatial mode \HG^o) in port B. 

In order to produce heralded single photons in arbi- 
trary spatial modes, one may apply a directional com- 
pression stress on the three-mode fiber with two plates 
as shown in Fig [TD1 (b). Such a compression breaks the 
cylindrical symmetry of the fiber medium which causes 
the non-cylindrically symmetric HG\q and HGqi modes 
to experience differing phase velocities as they propagate 
though the fiber, resulting in a nonzero relative phase 
at the fiber output that is controllable by the magni- 
tude and direction of the applied stress [21]. That is, the 
compressor acts analogously to a fractional wave plate for 
polarization, with principal axes aligned parallel and per- 
pendicular to the direction of compression, and with the 
degree of phase retardation proportional to the amount 
of stress. In the case illustrated, the biphoton mode 
function propagating in the fiber (before compression, 
but after being filtered) is \iPf) — \HG00) \HG^o) + 
I if G450 ) I H Goo) , the direction of compression is along the 
y axis, and the magnitude is such that so that the HG\q 
and HGqi mode components experience a relative phase 
factor of e % 2 = i. Therefore, the spatial part of the single- 
photon mode function |ifG45o) = |ifGio) + |ifGoi) be- 
comes ll/G^ 1 ) = \HGio)+i \HG01) after passing through 
the compressor, so that the heralded single-photon mode 
function at Sagnac output port B is both polarization 
entangled and in a well-defined first-order spatial mode: 

\i> B ) = \LG+ 1 )®(\H)\V) + \V)\H)) (13) 

We note that to cover the entire Poincare sphere of 
first-order transverse spatial modes [7], two such succes- 
sive compressors are needed with one applying twice the 
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amount of pressure as the other, in analogy with the re- 
quirement of a half- wave plate followed by a quarter- wave 
plate in order to turn linearly polarized light into light of 
an arbitrary polarization state. Also, the fiber compres- 
sor as shown in Fig. [10] (b) can be used in conjunction 
with the setup in Fig. \W\ (a) in order to produce, for 
example, LG-entangled bell states. Finally, we note that 
the above states can also be created by pumping the crys- 
tal with the appropriate higher order mode as opposed 
to using fiber compressors after the fact, since our Sorter 
can filter our any unwanted contributions to the bipho- 
ton mode function. However, in many cases, pumping 
the crystal with a Gaussian mode is preferable. 



CONCLUSIONS 
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We have demonstrated experimentally a stable out-of- 
plane Sagnac interferometer which sorts photonic trans- 
verse spatial modes according to two-dimensional par- 
ity. We performed 2-D parity sorting on Hermite- Gauss 
modes of definite parity and on the mixed-parity output 
field from a three- mode fiber. This interferometer is the 
first phase-stable device capable of measuring the trans- 
verse spatial mode parity of single-photon states, and 
is in principle 100% efficient. Because of this stability, 
many interferrometers can for the first time be practically 
cascaded in order to measure the orbital angular momen- 
tum of single photon states, without the use of holograms 
which are lossy in practice. The Sagnac can also be used 
to achieve spatial mode sorting of the output of a three- 
mode fiber into the fundamental Gaussian mode and the 
first-order mode family. This in turn makes feasible the 
generation of both HG-entangled Bell states and heralded 
single photons in arbitrary first-order spatial modes, cov- 
ering the entire Poincare sphere of first-order tranverse 
modes. In contrast to standard holographic techniques 
which path-separate spatial modes by converting them to 
the fundamental Gaussian, our Sagnac sorts modes while 
leaving their spatial mode functions unchanged, which 
may be desirable for certain applications. 

The 2-D parity sorting interferometer belongs to a 
larger class of interferometers which differ according to 
the relative rotation \I> of their constituent propagating 
and counter-propagating beams [2]. For example, the 2-D 
sorter (with \I/ = 180°) can distinguish between the zero 
order (ifGoo) and first order (ifGoi an d HG\q) modes, 
but not between the zero order and second order (HG02 5 
an d HG\i) modes. Conversely, a Sagnac with 
\I/ = 90° can distinguish between zero order and second 
order modes, but not between zero and first or third order 
modes. In general, the various interferometers discussed 
in this paper can be stably cascaded so as to sort photons 
belonging to any Poincare sphere order, which enables 
the manipulation of multi-dimensional qudits encoded in 
orbital angular momentum. 

Two-dimensional transverse spatial parity provides a 
novel type of qubit encoding that is based on continuous 



FIG. 11: A device that imparts a relative phase shift between 
counter-propagating fields, consisting of a tilt able birefringent 
waveplate (WP) with identical lengths of Faraday glass (FG) 
on each side, distinguished only by the direction of the ap- 
plied magnetic fields (B) permeating them, (a) A vertically 
polarized right-propagating photon becomes aligned with the 
fast axis of a birefringent waveplate (axis is denoted by the 
bold line), which induces a phase shift, (b) A vertically polar- 
ized left-propagating photon becomes aligned with the slow 
axis of the waveplate, thereby experiencing an unequal phase 
shift when compared to the forward- propagating beam. 

photonic degrees of freedom, in addition to the scheme 
based on one-dimensional parity previously realized by 
[6[. Unfortunately, to our knowledge it is not possible 
to design a Sagnac interferometer that sorts according 
to 1-D parity, since the relative transverse rotation upon 
which our Sagnac depends treats both spatial variables 
on an equal footing. Therefore, it appears that encodings 
based on 2-D parity can be made more robust than those 
based on 1-D parity due to the phase-stable manipula- 
tions of qubits enabled by this work. 
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APPENDIX A: MEASURING THE ODD- VALUED 
OAM OF SINGLE PHOTONS 

In this appendix, we propose an optical device that 
imparts a variable phase shift that is dependent upon 
whether the photon is propagating forwards or backwards 
through the device, as discussed in Section Hfl 

The device consists of a tiltable birefringent waveplate 
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surrounded on both sides by identical lengths of Faraday 
glass, distinguished only by the direction of the applied 
magnetic fields (B) permeating them as shown in FigfTTl 
The uniform magnetic fields cause the Faraday glass to 
act as a polarization rotator via the Faraday effect [221 ]. 
with the sense of rotation depending of the direction of 
the applied field. Upon entering the device from the left, 
a right-propagating vertically polarized photon under- 
goes a polarization rotation such that upon entering the 
waveplate the polarization is aligned with the waveplate's 
fast axis. The waveplate then imparts a phase shift to 
the photon before its polarization is returned to its orig- 
inal state, as shown in [TT] (a). Conversely, a vertically 
polarized left-propagating photon undergoes a transfor- 
mation which leaves its polarization aligned with the 
waveplate's slow axis. Therefore, the left-propagating 
photon experiences a different phase shift than the right- 
propagating photon before its polarization is returned to 
vertical, as shown in[TT](b). The waveplate can be tilted 
about its fast or slow axes in order to vary this relative 
phase shift, similarly to the standard method of tilting a 
thin glass slide. By employing one of these phase shift- 
ing devices inside our Sagnac interferometer, one may 
effect variable phase shifts between the propagating and 
counter-propagating beams without changing the polar- 
ization state of the light. We conclude that our Sagnac is 
capable of sorting odd OAM- valued photons from one an- 
other without the use of comparatively lossy holograms, 
and can therefore provide cascaded phase-stable sorting 
and measurement of single photons possessing arbitrary 
absolute OAM values, as discussed in Section ITTl 



APPENDIX B: EFFICIENT SORTING OF 
SINGLE PHOTON STATES 

In this appendix, we show that by replacing beam split- 
ter BS2 in Fig. [6] with a Faraday isolator and appropri- 
ately adjusting the Berek compensators, 100% sorting 
efficiency of single photon states can in principle be ob- 
tained, as stated in Section Hill 

A Faraday isolator (shown in Fig. [T2|) consists of a 
piece of Faraday glass with a polarizing beam splitter 
(PBS) on each side. The application of a uniform mag- 
netic field causes the Faraday glass to act as a polariza- 
tion rotator via the Faraday effect [22[. The first PBS 
(PBS1) is oriented at a positive 45° angle with respect 
to the vertical, such that an incoming photon polarized 
along this same direction will be completely transmitted. 
Conversely the second PBS (PBS2) is oriented vertically 
(at a 45° angle with respect to the first), in order to com- 
pletely transmit vertically polarized photons. Therefore, 
the polarization state of an incoming photon initially po- 
larized at 45° with respect to the vertical will remain un- 
changed as the photon is completely transmitted though 
PBS1, but will then be transformed to vertical polariza- 
tion so that the photon is also completely transmitted 
through PBS2, as shown in Fig. [12] (a). However, an ar- 




^ (No output) 

(No output) 



(b) 

Port A 



(No output) 




FIG. 12: A Faraday isolator, consisting of a piece of Faraday 
glass (FG) with a polarizing beam splitter on each side (PBS1 
and PBS2). (a) An input photon (propagating left to right) 
polarized at positive 45° angle with respect to the vertical is 
completely transmitted and becomes vertically polarized, (b) 
A back-propagating (from right to left) output photon from 
port A of the Sagnac with vertical polarization is completely 
transmitted through PBS2, but is subsequently completely 
deflected from PBS1. 



bitrarily polarized photon propagating through the iso- 
lator in the opposite direction will have its horizontally 
polarized component deflected by PBS2, while its verti- 
cally polarized component will be transmitted and sub- 
sequently rotated by the Faraday glass to a negative 45° 
angle with respect to the vertical, so that this remaining 
component is completely deflected by PBS1. 

Applying this to our current problem, one finds that 
although our backwards-propagating beam from output 
port A will be completely deflected as desired upon em- 
ploying a Faraday isolator, it would also in general be 
separated into two different paths by the two PBS's. Al- 
though one can easily recombine these beams using a half 
wave plate and an additional polarizing beam splitter, 
such a separation and recombination is undesirable as 
it introduces the same phase noise and drift that gives 
our setup an advantage over previous Mach-Zehnder- 
based 2-D sorter implementations. In our experiment, 
we can overcome this difficulty by adjusting the Berek 
compensators such that given a vertical input polariza- 
tion, the output polarizations of both port A and port 
B are vertical. Therefore, with the compensators ap- 
propriately adjusted, we conclude that back-propagating 
photons from port A will be completely transmitted by 
PBS2 and completely deflected by PBS1, as shown in 
Fig. [12] (b). In light of this, it is apparent that PBS2 of 
Fig. [12] is actually unnecessary, however we have included 
it because of the practical consideration that Faraday 
isolators are readily commercially available in the form 
discussed above. We conclude that through the use of 
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a Faraday isolator in conjunction with the Berek com- the input path, so that our 2-D parity sorting Sagnac is 
pensators, 100% of the output field from port A can be in principle 100% efficient with respect to both output 
deflected into a single beam path that is distinct from ports. 
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